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© In the encoding of a frame (P3) in a digitised 
motion picture sequence, the frame is divided into 
blocks (A) of pixels and a motion vector (VA) found 
to predicts the contents of each block A from a block 
(A') in a previous frame (10). Rather than evaluate 
many thousands of candidate vectors to find the 
best motion vector (VA), relatively few candidate 
vectors are evaluated, based on starting vector 
which is the estimated motion vector for a neigh- 
bouring block (LB(A)) in the same frame. To allow 
parallel processing of several pixel blocks, all the 
blocks of one row in the current frame are assigned 
to one processor (DSP), and a starting vector is 
obtained from a block (LB(A)) previously processed 
by the same processor. Further starting vectors are 
obtained if desired from neighbouring blocks (UR(B)- 
,DR(B)) in neighbouring frames of the sequence. The 
candidate vectors based on a starting vector may be 
very few in number, or may fill an area a few percent 
of the full search area. The method may favour 
candidate vectors which are more economical to 
encode. 
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The invention relates to a method of estimating 
motion vectors for blocks of pixels in a frame of a 
digitised motion picture sequence, the method 
comprising, for each given block of pixels perform- 
ing a vector search based around at least one 
starting vector derived from a previously estimated 
vector for a pixel block which is a neighbour in 
space and/or time of the given pixel block. 

The techniques of motion vector estimation can 
be used to improve compression of motion picture 
data by reducing the temporal redundancy in the 
picture information, in a manner analogous to that 
in which well-known differential coding and trans- 
form coding techniques eliminate spatial redun- 
dancy for the compression of still picture data. 

The International Standards Organisation has 
issued a Committee Draft Standard ISO 11172, 
well-known the art as the MPEG specification and 
incorporated herein by reference. MPEG specifies 
a format for the encoding of compressed motion 
picture data, including motion vectors. An overview 
of the MPEG coding standard is provided by D. Le 
Gall in 'MPEG: A Video Compression Standard for 
Multimedia Applications', Comm. of the ACM, Vol. 
34, No. 4, April 1991. MPEG encoding allows the 
reproduction of full-motion video from a data 
stream of only 1.2 Mbits per second, a rate low 
enough for the reproduction of video from MPEG 
data stored on a Compact Disc (CD-ROM), or re- 
ceived via some other data channel of similar ban- 
dwidth. 

The MPEG specification provides for motion 
compensation by dividing a frame into blocks of 16 
x 16 pixels, and assigning to each a motion vector 
in the range ±64 pixels in the x (horizontal) and y 
(vertical) dimensions, to half-pixel precision. This 
will be explained in more detail below with refer- 
ence to Figure 1 of the drawings. The MPEG 
specification does not specify how the motion vec- 
tors are to be obtained, only how they are to be 
represented. A full search for the most accurate 
motion vector in the range permitted would involve 
comparing the current pixel block with over 16000 
candidate pixel blocks, even at integer pixel preci- 
sion. Dedicated hardware can be provided for this 
task, but there is also a need for software im- 
plementations that can be used in advance of the 
general availability of such hardware, and which 
offer more flexibility to investigate improvements in 
the vector search and other aspects of the encod- 
ing strategy. 

A method as set forth in the opening paragraph 
has been described by G. de Haan and H. Muijgen 
in their paper 'Motion estimator for TV-picture en- 
hancement', presented at the Fourth International 
Workshop on HDTV, Turin, September 1991. That 
method has been developed for implementation in 
dedicated hardware and real-time operation, and 



involves the evaluation of relatively few candidate 
blocks. The method exploits the fact that, where 
object features moving in the motion picture se- 
quence are larger than a few pixels blocks, as is 

5 generally the case, there is a high degree of cor- 
relation between the motion vectors of neighbour- 
ing blocks. Thus, by careful selection of starting 
vectors from among vectors estimated for neigh- 
bouring blocks, and by the selection of candidate 

70 vectors predicted from the starting vectors, good 
quality motion estimation can be achieved with 
relatively few candidate block comparisons. While 
the known technique can be implemented in soft- 
ware to produce a quicker result than the full 

75 search, even using a high speed programmable 
digital signal processor chip (DSP) the encoding 
process is still inconveniently slow. 

A possibility for increasing the speed exists in 
theory in the use of available hardware systems 

20 comprising several DSPs operating in parallel, each 
processing a subset of the blocks in each frame. 
Unfortunately, this speed increase is not readily 
realisable with the method of de Haan and Huijgen 
because a given DSP processing a given block 

25 may have to wait while other DSPs process neigh- 
bouring blocks before the starting vectors for the 
given block are available. 

It is an object of the invention to enable im- 
proved motion vector estimation to be performed 

30 using programmable hardware, for example to al- 
low the use of parallel processing to achieve a 
faster result. 

The invention provides in a first aspect a meth- 
od as set forth in the opening paragraph, charac- 

35 terised in that the estimation of motion vectors for 
each frame is performed in parallel for at least two 
ordered subsets of the pixel blocks of the frame, 
the division and ordering of pixel blocks into the 
subsets and the selection of the neighbouring pixel 

40 block(s) being performed such that, where a neigh- 
bouring block whose estimated vector is to be used 
to predict that of the given block is selected from 
the same frame as the given block, the said neigh- 
bouring block is one of the same subset as the 

45 given block and is ordered in the subset so as to 
precede the given block. In this way, each subset 
can be processed at full speed, without waiting for 
results from the processing of blocks in the other 
subset(s). 

50 Each subset may comprise a row of pixel 
blocks, the processor estimating motion vectors for 
the pixel blocks of the row in a predetermined 
sequence, the motion vector for each given block 
(disregarding the first block in each row) being 

55 predicted by the processor from the motion vector 
estimated for a block which precedes the given 
block in the said sequence. The rows may be 
horizontal or vertical, but a horizontal bias to the 
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search may be advantageous since horizontal mo- 
tion is predominant in typical motion picture se- 
quences. A set of n processors may each take 
every n-th row of blocks in the image frame. Alter- 
natively, each processor may take the next un- 
processed row in the frame whenever it has fin- 
ished a row of pixel blocks. Many other arrange- 
ments are of course possible. 

In preferred embodiments of the method, for 
each given pixel block at least one further predic- 
tion of the motion vector is based on a motion 
vector estimated for a pixel block in a neighbouring 
frame in the motion picture sequence. By this 
means starting vectors can be chosen from neigh- 
bouring blocks effectively surrounding the given 
block, which improves the accuracy of the es- 
timated motion vectors when object features move 
over the given block from a wide range of direc- 
tions. 

It will be appreciated by those skilled in the art 
that the terms 'neighbouring' and 'preceding' as 
used herein with reference to frames, blocks or 
pixels cannot be interpreted as referring only to 
strictly adjacent frames, blocks or pixels. In an 
MPEG encoding embodiment, for example, it is 
clear that x neighbouring' frames may be separated 
at least by a number of frame intervals of the same 
order of magnitude as the batch size M. 

In a particular embodiment, the predictions for 
each given pixel block (disregarding the first block 
in each subset) are based on motion vectors es- 
timated previously for: 

- the said neighbouring block in the same sub- 
set as the given block; and 

- two blocks in a neighbouring frame occupy- 
ing spatial positions relative to the given 
block so as to define with the said preceding 
block in the same subset as the given block 
a triangular distribution generally surrounding 
the given block. 

In the prediction of a motion vector on the 
basis of the at least one starting vector, it may be 
for example that matching between the given block 
and a block in a previous frame is measured for an 
array of from 5 to 1 0 candidate vectors distributed 
around the starting vector. The candidate vectors 
may in particular be distributed within a few pixels 
of the starting vector, perhaps with a horizontal 
bias, and may include the starting vector itself. 

It will often be the case, in the known method 
as well as in a method embodying the invention as 
described above, that a good set of starting vectors 
is not available, and perhaps only a single default 
starting vector is available for the given block. Such 
cases arise at the edges of the frames of pixel 
data, and also when encoding the first predicted 
frame of a sequence, for example following a 
scene change. Although such cases are in the 



minority, a full vector, search may still be prohib- 
itively time consuming. 

The invention accordingly provides, in a sec- 
ond aspect, a method of estimating a motion vector 

5 for a given block of pixels in a frame of a digitised 
motion picture sequence, by performing an evalu- 
ation of a plurality of candidate vectors based ar- 
ound a starting vector, characterised in that the 
candidate vectors surround the starting vector and 

10 are evaluated in an ordered series of generally 
increasing distance from the starting vector. The 
ordered series of candidate vectors may form a 
spiral around the starting vector, for example. The 
area covered by the candidate vectors may be for 

?5 example an order of magnitude smaller than that 
required for a full search, but, if the vector es- 
timated for a neighbouring block is used as the 
starting vector, the process will tend to converge 
on an accurate motion vector within a few blocks, 

20 always assuming that the moving objects are larger 
than a few blocks. 

Moreover, the ordering of the candidate vectors 
ensures that the best candidate vector will tend to 
occur sooner rather than later in the sequence, and 

25 various checks can be implemented when evaluat- 
ing candidate vectors so as to terminate the evalu- 
ation of candidate vectors once the best candidate, 
or at least a "good enough" candidate has been 
found. 

30 It is a further object of the invention to improve' 
generally the motion picture quality encoded by a 
stream of motion-compensated image data. 

In a third aspect, the invention . provides a 
method of estimating a motion vector for a given 

35 block of pixels in a frame of a digitised motion 
picture sequence, which estimated vector is to be 
encoded in a data steam encoding the motion 
picture sequence, the method comprising: 

- generating a set of candidate vectors for 
40 evaluation; 

- evaluating the quality of prediction provided 
by each of the candidate vectors; and 

- comparing the results of said evaluations to 
select one of the candidate vectors as the 

45 estimated vector for encoding in a stream of 

data encoding the motion picture sequence, 
wherein the amount of data required to encode the 
estimated motion vector is dependent on the value 
of the estimated motion vector and wherein the 

so selection of the estimated vector from among the 
candidate vectors depends not only on the quality 
of prediction provided by each candidate, but also 
on the amount of data required to encode each 
candidate in a data stream, such that the candidate 

55 vector selected is not the one providing the highest 
quality of prediction, if it is determined that another 
one of the candidate vectors can be encoded in 
less data while providing a prediction quality which 
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is within a predetermined threshold of the highest 
prediction quality. 

This third aspect of the invention enables an 
overall improvement in the quality of motion pic- 
tures encoded within a bit stream at a limited data 
rate. If the quality threshold is set appropriately, the 
bit stream capacity which has been freed by more 
economical coding of the selected motion vector 
can be used to improve the encoding quality of 
other parts of the picture sequence, even if some 
of the freed capacity is required for increased error 
information resulting from the selection of a motion 
vector which does not quite give the best predic- 
tion for a given block. Of course, use of the inven- 
tion in its third aspect is not limited to any com- 
bination with the other aspects of the invention. 

As a particular example, the MPEG coding 
standard provides a special format for encoding 
blocks which have the zero vector as their motion 
vector. In some circumstances a zero vector does 
not have to be encoded at all. In certain picture 
sequences, including for example cartoon anima- 
tion sequences, it has been found that a significant 
amount of bit stream capacity is "wasted" when 
conventional encoders seek to encode "best" mo- 
tion vectors for a large number of blocks represent- 
ing, say, uniform blue sky. For such blocks, the 
only change from frame to frame in fact comprises 
only random noise superimposed upon the plain 
coloured area. It has been found that for these 
sequences the invention in its third aspect can lead 
to great savings in bit stream capacity, and hence 
great improvements in overall picture quality, when 
the quality threshold is set just high enough that 
the encoder will ignore the random noise motion 
information. 

The invention still further provides methods 
and apparatuses for performing motion vector es- 
timation, and for encoding and storing motion pic- 
ture sequences. These and further advantageous 
techniques and features will be apparent to the 
skilled reader in the following description of spe- 
cific embodiments of the invention. 

Embodiments of the invention will now be de- 
scribed, by way of example, with reference to the 
accompanying drawings, in which: 

Figure 1 illustrates the problem of motion vector 
estimation in the MPEG coding system; 
Figure 2 shows in block form a parallel DSP 
system suitable for implementing methods em- 
bodying the invention; 

Figure 3 illustrates the operation of the system 
of Figure 2 in accordance with the invention; 
Figure 4 illustrates the allocation of ordered sub- 
sets of pixel blocks between four DSPs in a 
method embodying the invention; 
Figure 5 illustrates the generation of a set of 
candidate vectors in accordance with a first 



search strategy; 

Figure 6 shows three alternative patterns for 

selection of starting vectors; and 

Figure 7 illustrates the generation of a set of 

5 candidate vectors from a starting vector in ac- 
cordance with a second search strategy. 
The embodiment to be described encodes 
digitised video images for storage and reproduction 
in accordance with the MPEG decoder specifica- 

w tion, cited above. Within this specification, there are 
many layers of compression encoding and much 
flexibility for improved encoding strategies to pro- 
vide greater image quality and/or compression, and 
the present description is concerned chiefly with 

75 improvements in the method of motion vector es- 
timation, which will not be defined in the standard. 
For full details of the MPEG coding system, refer- 
ence can be made to the MPEG specification or 
the article by Le Gall, both cited above. 

20 Figure 1 illustrates certain steps in the motion 
vector estimation for a sequence of five frames 
IO,B1,B2,P3,B4,P5 of pixel data which form part of 
a motion picture sequence to be encoded. MPEG 
defines three types of frames, with different de- 

25 grees of motion compensated compression. The 
sequence (any sequence) begins with an I- 
typeframe (intrapicture) IO, which is coded with 
reference only to itself, and relies on no information 
from neighbouring frames. The sequence continues 

30 with 'batches' of predicted frames, including B-type 
(bidirectional predicted) frames and P-type frames 
(predicted). 

The batch size M can vary typically from one 
to four frames, and each batch comprises (M-1) B 

35 frames followed by a P frame. Thus in the example 
sequence of Figure 1, frames B1.B2.P3 form a 
batch of M = 3 frames, while frames B4 and P5 
form a second batch with M = 2. It will be apparent 
that a batch with M = 1 comprises a single P-type 

40 frame. The arrows at the foot of Figure 1 indicate 
that each P frame is to be 'forward predicted' from 
the preceding I or P frame, while each B frame is 
to be both forward and backward predicted from 
the I and/or P frames on either side. 

45 The principle of motion compensated predic- 
tion will now be described briefly, with reference to 
the forward prediction of frame P3. Each frame to 
be predicted is divided into an array of blocks of 
16x16 pixels, represented (on an enlarged scale) 

50 by the block A in frame P3. Having encoded the 
starting or reference frame IO for the sequence, the 
block A is not directly encoded, but rather a motion 
vector VA is. encoded which predicts that the block 
A is the same as a 16x16 pixel block A' in the 

55 reference frame IO. Data compression is thereby 
achieved by temporal redundancy reduction. Since 
the two blocks A and A 1 will not generally be 
exactly the same, the difference (prediction error) 
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is then also coded with spatial redundancy reduc- 
tion by the discrete cosine transform technique 
(DCT), working on smaller blocks of 8x8 pixels. 

Thus the first problem faced by the MPEG 
coder is that of motion vector estimation for each 
block in the current frame P3: how to identify the 
block A 1 in the reference frame (10) which best 
matches the current block A, and hence identify 
the best motion vector VA? MEPG allows motion 
vectors to take a limited range in x (horizontal pixel 
dimension) and y (vertical) respectively where in 
the latest specification these ranges Rx and Ry are 
both ±64 pixels. This creates a limited forward 
search area FSA, shown in Figure 1, in which to 
locate the best matching block A\ but a full search 
even of this area would require the evaluation of 
many thousand candidate blocks. For each can- 
didate block, 16x16 = 256 pixel comparisons must 
be made, and the squared pixel differences accu- 
mulated (for example) to arrive at a mean squared 
difference value by which to measure the matching 
of the candidate block. While a full search may be 
possible with dedicated hardware, it is presently 
desired as mentioned in the introduction above to 
implement motion compensated image coding in 
software, using readily available processing hard- 
ware. 

Outside the field of MPEG encoding, de Haan 
and Huijgen, cited above, have developed motion 
vector estimation methods which measure match- 
ing for relatively few candidate blocks.based on the 
observations (i) that there is generally a correlation 
between motion vectors of neighbouring blocks and 
(ii) that subjective image quality can best be main- 
tained if the motion vector field is smooth, rather 
than accurate. In other words, there is no great loss 
in subjective picture quality if the estimated motion 
vectors are constrained to vary only a little from 
block to block within a frame. Thus, within the large 
search area FSA, only a few candidate vectors are 
used for matching, and these candidate vectors are 
based on starting vectors estimated previously for 
blocks neighbouring the current block, spatially 
and/or temporally. 

While the techniques of de Haan and Huijgen 
can be adapted to a software implementation, 
which operates many times faster than one with a 
full-search, this is still several hundred times real 
time, even with a high performance processor. The 
present inventors have looked to systems of par- 
allel digital signal processor chips (DSPs), such as 
the ANDROX ICS-400 image processor board, 
available from ANDROX in Boston, Massachusetts, 
and illustrated in block form in Figure 2. 

The system of Figure 2 comprises a circuit 
board housing four processors, referenced DSPO to 
DSP3, which can be controlled by a host computer 
system via a VME bus interface unit 200. The 
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board can be housed for example inside a Sun 
Microsystems or similar workstation-type computer, 
which provides a convenient and familiar develop- 
ment environment. The VME interface is connected 

5 to a high speed bus 202, through which flows all 
information processed by the DSP system. A bus 
controller 204 is provided to manage such oper- 
ations. A memory 206 is connected to the bus for 
the storage of the video data frames being en- 

70 coded. 

Each processor DSPO-3 is connected to the 
bus 202 by a pair of cache memories CMO and 
CM1. Each cache memory CMO or CM1 can in 
practice hold the current block of 16 x 16 pixel 

15 values, together with the pixel values of the search 
area FSA and other data necessary to estimate a 
motion vector for the block, although to this end 
the range (Rx.Ry) of candidate vectors is limited in 
the present embodiment to (±48,±48). In operation, 

20 one cache of the pair is connected to the DSP at a 
time, while the other cache can be unloaded and 
reloaded with data for processing via the bus 202. 
In this way, a constant supply of fresh data is 
available to keep the respective DSP working full- 

25 time. Each DSP further has a pair of program 
memories PMO and PM1 . Memory PMO is used to 
store the program for the DSP, while memory PM1 
is used to store temporary values and data which is 
common to several blocks. 

30 Figure 3 shows briefly in block form how the 
task of MPEG coding has been partitioned in the 
present embodiment. Original pixel data for a batch 
of M frames (plus reference frame) is received at 
IMIN, and a first module INTV performs the search 

35 for motion vectors to integer precision for the whole 
batch (for example frames B1.B2.P3), using exclu- 
sively the original pixel data. Only this first module 
INTV need be described in any detail for an under- 
standing of the present invention, but the remaining 

40 modules will be described briefly as follows. Four 
modules HLFV/MDEC/FDCT, Q/RLC, RLD/DEQ and 
IDCT/REC form a main loop of the encoding pro- 
cess. The module HLFV/MDEC/FDCT receives the 
integer motion vectors and calculates predicted 

45 blocks and the resulting prediction errors, making 
mode decisions and performing the forward DCT 
compression in accordance with MPEG published 
recommendations. In this process, the motion vec- 
tors are Tine tuned' to half-pixel precision. 

50 Module Q/RLC quantises and run-length codes 
the DCT information and passes these to a module 
BITS/REG which generates the final MPEG bit 
stream. The module BITS/REG also provides regu- 
latory feedback to the quantiser module to maintain 

55 the output bit rate at the desired value (typically 1.2 
Mbits per second). Returning to the main loop, 
modules RLD/DEQ and IDCT/REC reconstruct the 
pixel data actually encoded, which is fed back to 
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allow module HLFV/MDEC/FDCT to make correct 
mode decisions, locate the best half pixel vectors 
and so forth. Finally a module DIS allows display of 
the original or reconstructed pictures. 

In the implementation of these modules, the 5 
parallel system of four DSPs can only provide a 
four-fold increase in processing power over a sin- 
gle DSP if all four DSPs are busy all the time. To 
this end, the ANDROX system is interrupt driven so 
that each DSP operates independently, but this w 
creates problems for the implementation of the 
known methods of vector estimation in module 
INTV, because of the need to predict starting vec- 
tors from the vectors of neighbouring blocks. If 
neighbouring blocks are assigned to different 15 
DSPs, and particularly since the estimation may 
take a shorter time for some blocks than others, 
there will be times when DSP2, say, requires a 
starting vector which has yet to be calculated by 
DSP3. Because of this, much of the potential for 20 
parallel processing would remain unexploited. 

Figure 4 shows how the blocks of a frame can 
be partitioned into ordered subsets, as part of a 
solution to the above problem. In this embodiment, 
the blocks are divided into their horizontal rows, 25 
and all the blocks of a row are processed, in left to 
right order, by a single DSP. In the particular 
embodiment illustrated, every fourth row is pro- 
cessed by the same DSP, with the rows relevant to 
DSPO being highlighted in the Figure. In connection 30 
with the alternation of usage of the CMO and CM1 
caches, alternate blocks are shaded differently to 
illustrate which block is processed in which cache. 
Processing the blocks in horizontal rows is pre- 
ferred for most natural image sequences, where 35 
horizontal motion predominates. There is in princi- 
ple no reason why the blocks should not be pro- 
cessed in vertical rows (columns), diagonal rows or 
other sequences. 

It will now be explained how the module INTV 40 
estimates integer motion vectors, firstly for the 
blocks of the first batch of frames B1,B2,P3, repre- 
sented by block A in Figure 1, and then for the 
blocks of subsequent batches, represented by 
block B in frame P5. Two alternative vector search 45 
strategies are embodied in the implementation to 
be described, depending on whether the reference 
frame (I- or P-type) being used for the current 
frame has itself been motion predicted. (It may be 
noted here that I frames are often encoded in a so 
sequence to provide an 'entry points' into the se- 
quence, even though these frames can be motion 
predicted; for such frames it is in fact advanta- 
geous to find motion vectors to enable use of the 
second search strategy, even though these vectors 55 
will not be encoded in the MPEG bit stream). Both 
strategies rely on the choice of a good starting 
vector for the search, from which is generated a set 
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of candidate vectors, pointing to a corresponding 
set of candidate pixel blocks within the reference 
frame. 

Referring now to the vector search for block A, 
in accordance with the scheme of Figure 4, the row 
of blocks containing block A are all processed in 
left to right order by DSP3. Thus it is guaranteed 
that the block LB(A) immediately to the left of the 
current block A on the current frame P3 will have 
been processed already by DSP3, and an es- 
timated motion vector will thus be available which 
predicts the block LB(A) from some block in the 
reference frame IO. This vector can be used as 
one starting vector on which to base a search for 
the motion vector that will predict block A from a 
block A' in frame IO. On the assumption that no 
motion vectors have been found for the reference 
frame IO, nor yet the frames B1 and B2, a limited 
area-search strategy is used for the blocks of 
frame P3, as illustrated by Figure 5. 

Figure 5 shows a spiral search pattern which, 
on the integer pixel grid, extends ±15 pixels in the 
horizontal direction (Ax) and ±8 pixels vertically 
(Ay). This horizontal bias is provided in recognition 
of the fact that horizontal motion predominates in a 
typical motion picture sequence, as mentioned 
above. The program memory PM1 of the relevant 
DSP contains an ordered list of offset pairs (Ax.Ay) 
which define a set of candidate vectors '0' relative 
to a starting vector, indicated by ' + ' at offset (0,0). 
The sequence of candidate vectors 'o' spirals gen- 
erally outward from the starting vector at ' + 1 to fill 
the search area up to (Ax,Ay) = (±15,±8). In this 
way, even with a single starting vector, 527 can- 
didate vectors V are defined. Each candidate vec- 
tor in turn identifies a candidate block in the refer- 
ence frame IO to be compared with the current 
block A which it is desired to predict. 

For each block in the fourth row of blocks in 
the frame P3 to be predicted, the DSP must there- 
fore identify a starting vector. For the first block in 
a row, the zero vector (0,0) is used as a default 
starting vector, but thereafter the vector just es- 
timated for the block LB(A) to the left of the current 
block is used as the starting vector. This brings 
important benefits since in most cases there will be 
a high degree of coherence between neighbouring 
vectors in the vector field. Having identified the 
starting vector ' + \ the stored sequence of offsets 
can be used to generate the set of 527 candidate 
vectors V in the spiral sequence. 

For each candidate vector 'o\ the current block 
A in frame P3 is compared with a candidate block 
in reference frame IO whose spatial position (x.y) is 
given by the candidate vector relative to the spatial 
position (x,y) of the origin of the current block. The 
comparison is performed pixel by pixel in principle 
for all 256 pixels in the block, accumulating to a 
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mean squared difference (or other value measuring 
the degree of similarity between the blocks), until 
the best matching candidate block is identified. It 
may be noted, however, that a candidate block can 
in principle be rejected as soon as the accumu- 
lated difference exceeds the difference found for 
the best candidate block tried so far, resulting in a 
considerable saving of time. Since checking the 
accumulated squared difference against the best 
so far will generally occupy at least one cycle of 
the DSP's time, to make such a check after every 
pixel would slow down the entire matching process. 
A good compromise has been found in an embodi- 
ment where the accumulated squared difference is 
checked against the best so far only at the end of 
each row of 16 pixels in the block. 

While only a few per cent of the possible 
motion vector area FSA (Figure 1) is explored by 
this spiral search pattern, it is found that the es- 
timated vectors converge upon the true motion 
vector within a few blocks of the picture's edge. It 
may be noted that inaccuracies of the motion vec- 
tors do not directly reduce the quality of the cor- 
responding blocks in the image when decoded and 
displayed. Rather, inaccurate motion compensation 
merely reduces the degree of compression that 
can be achieved by the next stage in the MPEG 
coding process, namely spatial redundancy reduc- 
tion by DCT. 

Although far fewer than the number of can- 
didate blocks involved in a full search of the area 
FSA (Figure 1), 527 candidate blocks per block of 
the current frame is still a heavy computational 
task. The ordering of the candidates into a spiral 
brings potential savings, however, on the assump- 
tion that motion vectors are generally coherent 
from one block to the next. Because of this coher- 
ence, the starting vector chosen will, more often 
than not, be close to the desired motion vector for 
the current block. The ordering of candidates ac- 
cording to distance from the starting vector means 
that the best matching candidate block will tend to 
be found sooner rather than later in the sequence 
of 527 candidates. 

By incorporating in the DSP programming a set 
of test conditions designed to detect the best can- 
didate as soon as it appears, the search can be 
terminated, on average, after relatively few can- 
didate blocks have been compared to the current 
block. One such condition is the setting of a thresh- 
old vaJue of the accumulated mean squared dif- 
ference. Below this threshold (which allows, among 
other things, for noise in the pixel data), the blocks 
are judged a 'perfect' match, and the spiral search 
is terminated. Thus, for example, the generation 
and evaluation of candidate vectors 'o f might end 
with the discovery of a 'perfectly' matching can- 
didate block corresponding to the starting vector 



plus the offset (Ax.Ay) = (3,-2), as illustrated in Fig- 
ure 5. Up to this point, only 48 candidate blocks 
have been compared. 

Around the location of the best integer can- 

5 didate vector 'o\ represented at (Ax.Ay) = (3,-2) in 
Figure 5, nine crosses Y indicate the nine can- 
didate vectors that will be explored to half-pixel 
precision by the module HLFV/MDEC/FDCT, de- 
scribed above with reference to Figure 3. 

10 Clearly, the order of candidates can be modi- 
fied simply by modifying the list of offsets (Ax,Ay) 
stored in the DSP program memory PM1. The 
range of offsets Ax and Ay can be expanded or 
decreased as desired, and circular or elliptical spi- 

;s rals or other sequences can be implemented, offer- 
ing some improvement over the square and trun- 
cated spiral of Figure 5. 

Figures 6 and 7 illustrate the second vector 
search strategy, to be applied, by way of example 

20 to blocks in frame P5, represented by the block B. 
The strategy to be used for frame P5 compares far 
fewer candidate blocks even than the limited area 
search strategy described above with reference to 
Figure 5. This second strategy relies heavily on the 

25 choice of a few starting vectors, at least one of 
which should be a good prediction of the motion 
vector for the current block. Coupled with the rule 
that the vector field is better smooth than accurate, 
even these few candidates are found to allow a 

30 'good' yet extremely fast motion vector estimate. 

Turning to Figure 1, the block LB(B) imme- 
diately to the left of the current block B is again 
known to have been processed already, because of 
the ordered allocation of the whole row of pixel 

35 blocks to a single DSP. The vector just estimated 
for block LB(B) can therefore be used as one 
starting vector. Unfortunately, no other neighbour- 
ing blocks in the current frame P5 are known to 
have been processed at this time. The strategy 

40 implemented therefore looks to the neighbouring 
blocks which are known to have been processed, 
namely neighbouring blocks in neighbouring 
frames, such as the current reference frame, the 
previous P frame P3. 

45 For the representative block B, two blocks UR- 
(B) and DR(B) are shaded in frame P3, these being 
spatially diagonally up-and-right and down-and- 
right neighbours of the block B respectively, as 
well as temporal neighbours of the current block B 

so (frame P5). The two neighbouring blocks provide a 
further two starting vectors for use in estimating the 
motion vector of block B. 

It may be noted at this point that such starting 
vectors may need to be scaled versions of the 

55 vectors estimated for the neighbouring blocks. For 
example the vector VA predicting block A in frame 
P3, say, from block A' in frame 10 represents 
motion of some image feature over a time period of 
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three frame intervals (M = 3), whereas a vector pre- 
dicting block B in frame P5 from a block B' (not 
shown) in frame P3 represents motion over two 
frame intervals only. Thus the vector estimated 
previously for block* UR(B) must be scaled by a 
factor of 2/3 to obtain a good starting vector for the 
vector estimation of block B, and so forth. (If back- 
ward prediction vectors, found for blocks in a B 
frame, are being estimated or used as starting 
vectors, negative scaling factors will also be re- 
quired). 

Figure 6 A shows the pattern formed by the 
chosen neighbouring blocks around the current 
block, with different hatching to distinguish blocks 
which are spatial neighbours in the same frame 
(t = T) and blocks which are both spatial and tem- 
poral neighbours (t*T). The chosen neighbour 
blocks generally form, a triangle around the current 
block, with a slight horizontal bias to the distribu- 
tion. The horizontal bias is in recognition of the fact 
that in most natural motion picture sequences, hori- 
zontal motion predominates. Choosing neighbour 
blocks to the left and right of the current block 
helpsto ensure that at least one of the starting 
vectors will be suitable when an object moving into 
the current block comes from the left or the right. 

Figure 6B shows an alternative pattern of 
neighbouring blocks to yield three starting vectors, 
and Figure 6C shows a pattern yielding four start- 
ing vectors, without horizontal bias. Always in Fig- 
ure 6 the block LB to the left of the current block is 
in the current frame (t = T). Of course other patterns 
can be readily envisaged by those skilled in the art, 
each having its own merits and demerits, depend- 
ing on the motion picture content, and the pattern 
can be changed by changing the contents of the 
program memories PMO, PM1 for each DSP. The 
pattern of Figure 6A has been found most satisfac- 
tory in the present embodiment, for a range of 
picture sequences. Of course, yet different patterns 
will be required where the allocation of blocks to 
the DSPs is different to that of Figure 4. The block 
rows might run vertically or even diagonally across 
the picture, or the subsets might be convoluted 
and/or interleaved for whatever reason; the princi- 
ples disclosed herein can still be applied. 

Having obtained three starting vectors for the 
current block A, the DSP in the present embodi- 
ment evaluates the block matching for a group of 
up to 21 candidate vectors based upon these three 
starting vectors. 

Figure 7 shows how seven candidate vectors, 
marked by circles 'o\ are defined by a predeter- 
mined pattern of offsets (Ax,Ay) relative to the 
starting vector, marked ' + \ The pattern of seven 
candidates 'o' includes the starting vector 1 + 1 itself, 
and again has a horizontal bias in recognition of the 
predominance of horizontal motion in the human 



environment. By applying this pattern of seven 
offsets to each of the three starting vectors, a set 
of 21 candidate vectors is obtained. Before any 
time-consuming block matching calculations are 
5 performed, however, the candidate set is checked 
to eliminate duplicates and candidates outside the 
legal range. 

Once the best matching candidate vector 'o' 
has been found, again a final comparison is per- 

w formed in module HLFV/MDEC/FDCT (Fig. 3) with 
nine half-pixel candidate vectors, marked as V 
around the rightmost candidate vector V in Figure 
7. This leads to a final encoded motion vector for 
the current block, estimated to half-pixel precision 

75 and derived using reconstructed pixel data. 

It will be apparent to the reader that evaluating 
at most 21 candidate vectors instead of the thou- 
sands required for a full search of the area FSA 
(Figure 1) achieves a significant time saving, but 

20 with the disadvantage that the vector estimated for 
the current block may not in fact be the most 
accurate one. The subjective quality of images 
decoded after encoding by this method is good, 
however, because of the above mentioned observa- 

25 tion by de Haan and Huijger that the motion vector 
filed should be above all smoothly-varying, rather 
than strictly accurate. Moreover, when object fea- 
tures moving in an image are many pixel blocks in 
size, the slowly varying vector field will tend to 

30 converge upon the most accurate vector within a 
few blocks after entering the object. 

Figure 1 also shows how, on the assumption 
that motion vectors have been estimated already 
for the interpolated frame B4, as well as the refer- 

35 ence frame P3, a further four possible starting 
vectors are obtainable from by scaling the forward 
and backward vectors estimated for the UR and DR 
neighbouring blocks in frame B4, shaded with dot- 
ted lines in Figure 1. Thus a total of up to six 

40 starting vectors are available in addition to the 
vector estimated for the preceding block LB(B). A 
further four or eight starting vectors could be avail- 
able for the last frame in a batch of three or four 
frames. Rather than evaluate a set of seven can- 

45 didate vectors for every one of these starting vec- 
tors, the present embodiment in fact first evaluates 
only the starting vectors themselves, before for- 
ming candidate sets around the starting vectors 
which give the best matching candidate blocks. 

so Thus, for example, candidate sets according to the 
pattern of Figure 7 may be formed on the basis of 
vectors estimated for the block LB(B), the block 
UR(B) from frame P3 and the backward prediction 
vector for the DR block in frame B4, all scaled as 

55 necessary to account for the different time intervals 
spanned. 

In practice, the vector search for each batch of 
M frames is performed in the order indicated by 
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circled numbers by the curved arrows at the foot of 
Figure 1. By adopting this order, and with appro- 
priate scaling, forward and backward prediction 
vectors of neighbouring blocks in the interpolated 
frames B1 and B2 are all available as candidate 
starting vectors when estimating motion vectors for 
frame P3. Similarly, the forward vectors for frame 
B1 can be scaled and used as starting vectors in 
the forward vector search for frame B2. 

The remainder of the MPEG encoding process 
involves variation of many parameters such as the 
batch size M and quantisation scaling, in order to 
achieve a desired compromise between picture 
quality and data compression, as described in the 
MPEG specification and the article by Le Gall. 
These aspects are outside the scope and purpose 
of the present description, but it has been found 
that the complete MPEG coding of motion picture 
sequences including the motion estimation meth- 
ods described above can be achieved at only 90 
times real time using the four DSPs of the AN- 
DROX system, mounted in a Sun workstation. Thus 
a few minutes of full motion video can be encoded 
for storage on CD-ROM or similar on an overnight 
basis, using readily available hardware in a small, 
compact system. 

A further modification of the above strategies 
can be made when the encoding scheme results in 
some vector values being more economical to en- 
code than others. In the particular examples of the 
MPEG encoding format, motion vectors are en- 
coded in the MPEG bit stream by variable length 
codes, such that shorter vectors require fewer bits. 
Thus, when two candidate vectors have almost 
equal matching quality, the modified strategy will 
favour the vector which requires fewer bits to en- 
code. This may necessitate a little more error in- 
formation being encoded for a given block, but 
there will often be a net saving in encoded data 
which allows the image to be encoded with an 
overall improvement in quality. 

A particular feature of the MPEG variable 
length codes for motion vectors is that zero dis- 
placements Ax and Ay can be encoded in only one 
bit each. Moreover, at least in a p-type frame a 
block having zero motion (Ax=Ay = 0) can be cod- 
ed with no motion information at all. This is sig- 
nalled by the appropriate variable length code for 
"macroblock type" in the MPEG bit stream. Thus 
there is a great economy of data in encoding 
blocks which have zero motion from one frame to 
the next. Thus an overall improvement can be 
obtained when the evaluation of candidate vectors 
is modified to favour the zero vector, when the 
matching quality obtained is only slightly below 
that of the best candidate. 

Of course, in some sequences, the zero vector 
will rarely be among the best candidates. It may be 
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noted however, that certain types of motion picture 
sequence, for example animated cartoons, contain 
a large number of blocks with substantially zero 
motion. Moreover, large areas of these block com- 

5 prise uniform areas of colour, for example blue sky. 
Conventional motion compensation techniques will 
nevertheless estimate and encode motion vectors 
for these blocks, if only because they contain noise 
information which "moves" from frame to frame. 

w For encoding these sequences, it has been found 
that very significant quality improvements are af- 
forded by use of the invention in its third aspect: if 
the noise level is low, the zero vector will be just 
one of a large number of candidate vectors for 

75 which the matching quality is the same within a 
relatively narrow threshold. 

From reading the present disclosure, many 
modifications will be apparent to persons skilled in 
the art, in addition to those mentioned above. Such 

20 modifications may involve other features which are 
already known in the fields of picture encoding, 
motion estimation and digital signal processing and 
which may be used instead of or in addition to 
features already described herein. Although claims 

25 have been formulated in this application to particu- 
lar combinations of features, it should be under- 
stood that the scope of the disclosure of the 
present application also includes any novel feature 
or any novel combination of features disclosed 

30 herein either explicitly or implicitly or any 
generalisation thereof, whether or not it relates to 
the same invention as presently claimed in any 
claim and whether or not it mitigates any or all of 
the same technical problems as does the present 

35 invention. The applicants hereby give notice that 
new claims may be formulated to such features 
and/or combinations of such features during the 
prosecution of the present application or of any 
further application derived therefrom. 

40 

Claims 

1. A method of estimating motion vectors for 
blocks of pixels in a frame of a digitised mo- 

45 tion picture sequence, the method comprising, 
for each given block of pixels performing a 
vector search based around at least one start- 
ing vector derived from a previously estimated 
vector for a pixel block which is a neighbour in 

so space and/or time of the given pixel block, 
characterised in that the estimation of motion 
vectors for each frame is performed in parallel 
for at least two ordered subsets of the pixel 
blocks of the frame, the division and ordering 

55 of pixel blocks into the subsets and the selec- 
tion of the neighbouring pixel block(s) being 
performed such that, where a neighbouring 
block whose estimated vector is to be used to 
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predict that of the given block is selected from 
the same frame as the given block, the said 
neighbouring block is one of the same subset 
as the given block and is ordered in the subset 
so as to precede the given block. 5 

A method as claimed in Claim 1 wherein each 
subset comprises a row of pixel blocks, the 
processor estimating motion vectors for the 
pixel blocks of the row in a predetermined w 
sequence, the motion vector for each given 
block (disregarding the first block in each row) 
being predicted by the processor from the 
motion vector estimated for a block which pre- 
cedes the given block in the said sequence. 75 

A method as claimed in Claim 1 or 2 wherein 
for each given pixel block at least one further 
prediction of the motion vector is based on a 
motion vector estimated for a pixel block in a 20 
neighbouring frame in the motion picture se- 
quence. 

A method as claimed in Claim 3 wherein for 
each given pixel block in a row (disregarding 25 
the first block in each row) predictions are 
based upon the previously estimated motion 
vectors of the said preceding block in the 
same row and a plurality of blocks in the said 
neighbouring frame. 30 

A method as claimed in any of Claims 1 to 4 
wherein the predictions for each given pixel 
block (disregarding the first block in each sub- 
set) are based on motion vectors estimated 35 
previously for: 

- the said preceding block in the same 
subset as the given block; and 

- two blocks in a neighbouring frame oc- 
cupying spatial positions relative to the 40 
given block so as to define with the said 
preceding block in the same subset as 

the given block a triangular distribution 
generally surrounding the given block. 

45 

A method as claimed in any preceding claim 
wherein on the basis of the at least one start- 
ing vector, matching between the given block 
and a block in a previous frame is measured 
for an array of from 5 to 10 candidate vectors 50 
distributed around the starting vector. 

A method as claimed in any of Claims 1 to 5 
wherein on the basis of the at least one start- 
ing vector matching between the given block 55 
and a block in a previous frame is measured 
sequentially for an ordered series of candidate 
vectors surrounding the starting vector and 



with generally increasing distance from the 
starting vector. 

8. A method as claimed in Claim 7 wherein the 
evaluation of candidate vectors is terminated if 
a vector is found giving a matching better than 
a predetermined threshold. 

9. A method of estimating a motion vector for a 
given block of pixels in a frame of a digitised 
motion picture sequence by performing an 
evaluation of a plurality of candidate vectors 
based around a starting vector, characterised 
in that the candidate vectors surround the 
starting vector and are evaluated in an ordered 
series of generally increasing distance from 
the starting vector. 

10. A method as claimed in Claim 9 wherein the 
evaluation of candidate vectors is terminated if 
a vector is found giving a matching better than 
a predetermined threshold. 

11. A method as claimed in any of Claims 1 to 10 
wherein the selection of each estimated vector 
depends not only on the measured prediction 
quality (matching) of each candidate vector, 
but also on the relative economy with which 
the candidate vector can be encoded in a 
resultant data stream encoding the motion pic- 
ture sequence. 

12. A method as claimed in Claim 10 wherein the 
zero vector is most economically encoded and 
wherein a candidate vector which is zero will 
be favoured in the selection of each estimated 
motion vectors. 

13. A method of estimating a motion vector for a 
given block of pixels in a frame of a digitised 
motion picture sequence, which estimated vec- 
tor is to be encoded in a data stream encoding 
the motion picture sequence, the method com- 
prising: 

- generating a set of candidate vectors for 
evaluation; 

- evaluating the quality of prediction pro- 
vided by each of the candidate vectors; 
and 

- comparing the results of said evaluations 
to select one of the candidate vectors as 
the estimated vector for encoding in a 
stream of data encoding the motion pic- 
ture sequence, 

wherein the amount of data required to encode 
the estimated motion vector is dependent on 
the value of the estimated motion vector and 
wherein the selection of the estimated vector 
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from among the candidate vectors depends 
not only on the quality of prediction provided 
by each candidate, but also on the amount of 
data required to encode each candidate in the 
data stream, such that the candidate vector 5 
selected is not the one providing the highest 
quality of prediction, if it is determined that 
another one of the candidate vectors can been 
coded in less data while providing a prediction 
quality which is within a pre-determined 10 
threshold of the highest prediction quality. 

14. A method as claimed in any preceding claim 
wherein the motion vectors and prediction er- 
ror information are encoded in accordance with ts 
the MPEG standard as defined by the Interna- 
tional Standards Organisation. 

15. An apparatus for encoding a digitised motion 
picture sequence comprising: 20 

- means for dividing at least one given 
frame of the sequence into a plurality of 
blocks of pixels; 

- means for estimating by the method of 

any preceding claim a motion vector pre- 25 
dieting each block of pixels in the given 
frame from a block in a neighbouring 
frame and encoding said motion vector; 
and 

- means for generating and encoding pre- 30 
diction error information to for use with 

the motion vectors in reconstructing the 
encoded image. 

16. A method of storing a digitised motion picture 35 
sequence on a record carrier, comprising: 

• using an apparatus as claimed in Claim 
15 to encode the motion picture se- 
quence; and 

- storing the resultant encoded motion 40 
vector information and prediction error 
information on the record carrier. 

17. A record carrier wherein a motion picture se- 
quence has been stored in accordance with a 45 
method as claimed in Claim 16. 
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